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I. INTRODUCTION 



The theory of signal detection in additive white Gaussian noise (WGN) is well 
established. Optimum structures and their performance can be found in many 
textbooks [Ref 1; Chap. 4). Solution of the so-called Binary Hypothesis Testing 
problem allows one to obtain an optimum structure (receiver) that is capable of 
deciding with minimum probability of error, which of two possible signals, Sj(t) or 
s,^(t), was transmitted over an additive white Gaussian noise (AWGX) channel. Such a 
receiver is shown in Figure 1.1 




Figure 1.1 Structure of the Optimum Receiver, 
and its performance (i.e., probability of error, P^) is given by [Ref 2: Chap. 6], 



P 



e 



Po r f 'i' 

— ^ erfc [ 

- y2NoE(l-p) 



-b 




erf [ 



Y - E(1 -p) I 
72 No E(1 - p) 



( 1 . 1 ) 



9 



\\here 



pj = Probability that Sj(t) was transmitted, i =0,1 

Nq 2 = Power spectral density (PSD) level of the additive WGN 

E = J [ Sgit) + Sj"(t) 1 dt 

- T 

T = Length of the observation interval over which signals are received 
1 r 

P = j So(t)s,(t) dt 

t j 

Y = Threshold = — ^ Cn Xf,. 



( 1 . 2 ) 



(1.3) 



(1.4) 



The erfc(x) and erf(x) functions are defined by 
00 1 

erfc(x) = I e'b ^ d^ 

X 



1 ^ 2 , 

erf (x) = j e''’ " d^ 

-CO V 2Jt 



and Xq depends on pQ and pj primarily, but may also depend on costs associated with 
each decision (whether correct or incorrect). If these costs are included but are all 
assumed to be equal, and Pq = pj. then Xg = 1. 

The design of the receiver shown in Figure 1.1 does not take into account the 
possibility that in the transmission channel, in addition to the AWGX, another source 
of interference (perhaps intentional, such as a januner) may be present, with 
characteristics far difTerent from the assumed AWGN interference. 

In this thesis, we analyze the performance of the receiver of Figure 1.1, when in 
addition to the AWGN, an additive nonstationary Gaussian noise is present in the 
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channel. In Chapter II. this nonstationary Gaussian noise component is described, and 
its possible generation by an intentional interl'erer is analyzed. In Chapter HI. the 
performance of the receiver of Figure 1.1 is analyzed taking into account the presence 
of the nonstationary interference. It is demonstrated that the resulting performance is 
in all cases worse (i.e.. higher P^) than that predicted by Equation 1.1. (This is to be 
e.xpected. since the model that gave rise to Equation 1.1 is not optimum for any 
interference other than the AWGN). Chapter IV presents a modification to the receiver 
of Figure 1.1 that results in a structure that is optimum for the detection of binan.' 
signals in AWGN and additive nonstationary interference of the type modeled in 
Chapter 11. The performance of this modified structure is evaluated and compared with 
the results of Chapter III. In Chapter V several examples are presented for two 
important types of signaling schemes, namely Phase Reversal Keying (PRK) and 
Frequency Shift Keying (FSK) under various assumptions about the characteristics of 
the interference. The derived mathematical results are used to evaluate receiver 
performance in terms of probability of error. P^, as a function of signal to noise ratio, 
and interference to signal ratio. The results are then interpreted and conclusions are 
drawn. 



II. NOiNSTATIONARY INTERFERENCE MODEL 



As pointed out in Chapter I. the intentional interference nj(t) is modeled as a 

nonstationary Gaussian process having autocorrelation function R_ (t,T). A possible 

"... 

method for generating such a process is depicted in Figure 2.1, where a Gaussian 
process n„(t) is multiplied by the deterministic function q(t) to produce the process 
n.(t). 




Figure 2.1 Structure for the Generation of m(t). 

That is 

n.(t) = q(t)ng(t) 
so that 



R (t,t) = E{ n.(t)n.(T) ) = q(t)q(t)R (t,T) 



where 



R (t.t) = E{ n„(t)n (t) ). 

^ ^ a CO 

Note that m(t) is a Gaussian process since n„(t) is assumed to be Gaussian. If 
n (t) is white with unit power spectral density level, then R_ (t,T) = 6(t-T) so that 

c ^ U 

R (t.T) = q(t)q(T)6(t-T) = q-(t)6(t-t) = Q(t)6(t-T) (2.1) 

J 
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where 



Q(t) ^ q-(t) > 0. 
irrig(t) is a nonwhite process, then 

E(nj(t)nj(t+T)J = R^, (t.t + r)= q(t)q(t+ (r). (2.2) 

Implicit in these expressions is the assumed wide sense stationarity of iig(t). 
Observe from Equation 2.1 and Equation 2.2, that a(t) is a white (colored) 
nonstationary’ Gaussian process when n^(t) is a white (colored) Gaussian process. 

Since 



< R (t.t+t)> = (T) = 

J J 



Eim 

T-» ao 



1 



2T 



T 

I q(t)q(t + t)Rj^^ (r) dt 
-T ' 



= Q3(t)R 



n. 



(T) 



where 



Qa(^) = J (2.3) 

1 2 . 1 

-T 

the average power of n.(t) is given by 

Pnj ^ -^nj (0» = Qa(»»Rn„ (»)' (--» 

This expression has physical meaning only when R^ (0) is finite (i.e..n (t) is non- 
white). Note that when n^(t) is white (with unit PSD level). Qg(0) becomes the average 
PSD level of m(t). 

For most the work in this thesis, it will be assumed that the interference m(t) is a 
white nonstationary Gaussian process, with autocorrelation given by Equation 2.1. 
Several forms of q(t) used in the examples worked out in Chapter V are shown in 
Figure 2.2. 
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Figure 2.2 Several Examples of q(t) as a Function of Normalized Time. 
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III. RECEIVER PERFORMANCE ANALYSIS 



We now evaluate the performance of the receiver depicted in Figure l.I when 



I 



Sj(t) + njt) + m(t) w.p. pj 



r(t) = 



t e T, 



(3.1) 



So(t) + Po 



where i\.(0 represents the AWGN with PSD level Nq 2. and m(t) is the white 
nonstationary Gaussian noise with autocorrelation function given by Equation 2.1. 
Observe from Equation 3.1 that r(t) is a Gaussian random process, and its processing 
by the receiver of Figure 1.1 will result in G being a Gaussian random variable. 
Observe from Figure 1.1 that the receiver performs the test 



and if G > Y presence of Sj(t) is declared, whereas if G < Y- tlie presence of S(^(t) 
is declared. Whether Sj(t) or SQ(t) was transinitted, G is a Gaussian random variable 
having conditional means 



G = J r(t)s/t) dt + — J [ SQ-(t) - s,-(t) ] < Y 

T " T 



mj = E{G Sj(t) was transmitted] = — J [ Sj(t) - Sg(t) ]“ dt 

T 



(3.2) 



rrig = E{G;Sg{t) was transmitted] = - — j [ s^(t) - Sg(t) dt = - m^ (3.3) 

T 



and conditional variances 



var{G Sj(t) was transmitted] = E { ( J [ n^^.(0 + rij(t) | s^(t) dt |" ] 

T 
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where we have assumed that tuid m(t) are uncorrelated ^ero mean random 

process. Let 

*G- ~ probability density function ofG conditioned on Sj(t) being 

transmitted, i = 0.1 



so that 



Pg j = Pfdecide Sj(t) transmitted; SQ(t) was transmitted} Pq 
+ P(decide SQ(t) transmitted Sj(t) was transmitted} pj 

oo Y 

== Po J ^Go(?o)^?o + Pi 1 

y -'X) 

(the subscript s denotes that receiver is suboptimum for the noise model used), and 
from Equations (3.2) - (3.4) 



Po P, 

Pe,s - 7 — erfc [ (y - m^) V 2 (Tq ] + ^— erf [ (7 - nij); V 2 Oq ]. (3.5) 

Comparison of Equation 1.1 with Equation 3.5 is best achieved if we assume equal 
decision costs and equal prior probabilities pg and Pj, so that y = 0 (Under these 
assumptions Xg = 1 which results in y becoming 0). Thus from Equations (1.1) - (1.4) 



Pg = erfcl^Ed -p) 2Ngl (3.6) 

(note that E(1 - p ) ^ 0 ) and from Equation 3.5 

^e,s = erfc ( nij 7T <Tq ). 
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From Equation 3.2 and Equation 3.4 we obtain 



^ = [ -^ J Sj-(t) dt ]‘ [ ^ J Sj‘(t) dt + J Q(t)Sj‘(t) dt 

2 T ^ T T 



E(1 -P) 



J Q(t)Sj-(t)dt 



I - 



NqE(1 - p) + J Q(t)s^"(t) dt 
T 



and therefore 



Pes = erfc[yE(l-pXl-C)2N, 



(3.7) 



where 



C = j Q(t)Sj-(t) dt / [ XqE( 1 - P) + J Q(t)s/(t) dt ]. (3.8) 

T T 

Observe that C e [0.1] and therefore is the factor tltat causes an increase in probability 
of error due to the presence of m(t). If m(t) were not present, then Q(t) = 0 -♦ C = 
0 and Equation 3.7 would become identical to Equation 3.6. Thus we can study the 
performance degradation of the receiver by either evaluating or by evaluating C. 
as its size will dictate the increase in P^ Note that if Q(t) -» co, then C -> 1 and P^ ^ 
-♦ 1 2 as expected. It is worthwhile noting that if the nonstationary interference 
become (white) stationary so as to simply have the elTect of raising the AWGN level, 
then, with q(t) = K (a constant), from Equation 3.8 



C = 



K' J Sj-(t) dt 

T 






N 



^js,'(t)dt + K^Js,2(t)dt 
T T 



1 + K- 



(3.9) 
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Since the interference source will typically be an intentional jammer (in this case the 
jammer lacks sophistication), we can think of the K" (Nq 2 ) ratio as 



K' 



K-E 



Xo2 



){ 




(signal to noise ratio)(jammer to signal ratio) 



and thus 



P 



e.s 



erfc [ 



E(l-p) I 
2\ 1 + JSR-SNR 



(3.10) 



and 



C = JSR-SXR (1 + JSR-SXR) 



where 



SXR = E (Xo 2) 

JSR = K-E. 

We see that a large JSR value is required for the system to become useless as a 
receiver (namely Pg 5 12 ). It therefore appears that a jammer can use its available 

power more elTiciently by not spreading its power over large bandwidths. 
Consequently, jammer waveforms other than broadband noise power will be 
investigated in Chapter V. Xote furthermore that an increasing value of Xq in order to 
prevent JSR'SXR from getting too large is self defeating because as revealed by 
Equation 3.10. with JSR constant, 1 2 as Xq -* 

The derived equation for P^ ^ will be used in the sequel to analyze specific 
modulation schemes and choices of nonstationaiy waveform q(t). It will be 
demonstrated that by proper choice of q(t), the receiver of Figure 1.1 can be rendered 
inelTective without the use of a great deal of janmier power. 
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IV. RECEIVER PERFORMANCE ANALYSIS WITH INTERFERENCE 

STATIONARIZATION 



The optimum receiver for processing r(t) (with r(t) as given in Equation 3.1) is 
easily determined by realizing that the nonstationary interference can be stationarized 
by use of the (time-var\’ing) system depicted in Figure 4.1 




r(t) = Sj(t) + n^.(t) + nj(t) i = 0,1 

where n^,_.(t) and m(t) are described in the discussion following Equation 3.1, the output 
of the linear system depicted in Figure 4.1 is 

r'(t) = s/(t) + n'(t) i = 0,1 

where 



S;'(t) 



S;(t) 




n'(t) 



^ + n,(t) 

J ^ + q=(.) 



i = 0,1 



(4.1) 
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Evaluatine now the autocorrelation function of n't t). we obtain 



E[n'(t)n'(r),l = E( 



[njt) + nj(t)j [ n,Xt) + a(T) ] 



N,- 



V 






-W 



V 



+ q-(t) 



(t.t) + (t.T 



"w 




+ q-(t) 




q“(^) 



(4.2) 



where the uncorrelatedness of n^^.(t) ^ri<^ r^j(E) has been used in Equation 4.2. From 
Equation 2.1. we obtain 

Nn 

6(t - T) + q-(t)5(t - T) 

Rn' : = 5(t-t) 

yi^ + qn.) yy^ + q-(0 

and clearly, the output of the linear time var>ing system consists of a known signal 
s.'(t). i = 0,l, in AWGX of unit PSD level. The optimum receiver for deciding whether 
Sj'(t) or S(^'(t) was received in AWGX of unit PSD level is given by the receiver of 
Figure 1.1, with r(t) replaced by r'(t), Sj(t) replaced by 

Sj '( t ) = S ('( t ) - Sq '(0 

the bias term given by 



1 



•) 



and 



j' ( So'-(t) - Sj'-(t) ] dt 

T 



Y = Cn >-o. 
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The structure of this optimum receiver is shown in Figure 4.2. 













r(t\ ► 


1 


IIILsa , 


O < 






I L ) ^ 


- +q“(t) 




,-VTf^ 

bias . 

r 

Y 


comp, decision 




Sd'(t) 






^ t) = t) ~ = 




i = 0.I 






J 


+ q-(.) 



Figure 4.2 Structure of the Optimum Receiver with Nonstationar\’ Interference. 
Its performance is given by Equation l.I with 



Y = Cn Xf 



E replaced by E', where 



E' = -^ J'[So'“(t) + Sj'“(t)]dt = 



1 SQ-(t) + s^-(t) 



T 



T 



NV 



q“(t) 



(4.3) 



p replaced by p', where 



P' ^ j' I So'(t) Sj’(t)] dt = 



1 j dt 



N. 



+ q-(t) 



(4.4) 
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and Ng is replaced by 2 (since N g 2 = 1 in this case). Assuming equal prior 

probabilities and decision costs, we can obtain from Equation 3.6. with replacements 
indicated above, the performance of this optimum receiver, namely. 



Pg 0 = erfc [ VE'(1 - p') 4 ] 



(4.5) 



(the subscript o denotes that receiver is optimum for the model used). Since 



-p> ' t , 



N. 



0 T 



and from Equation 4.3 and Equation 4.4 



E'(l-P') 1 






^^0 T 1 + q-(tf^ 



dt 



(4.6) 



it is clear that 



E'(l-p')2 < E(l-p)Ng 

so that Pg Q > Pg. This means that the optimum receiver designed to operate in a 
nonstationary interference environmemt. performs worse that the optimum receiver 
designed to operate in the presence of WGN only, even though the former 
stationarized the interference prior to performing the (standard) correlation operation. 
This is not surprising because the former receiver is operating at a higher level of 
interference than the latter receiver. (Observe that with q(t) = 0. Pg = Pg q 
e.xpected). Furthermore, we must have o ~ ^e s' demonstrated by use 

of the Cauchy-Schwarz inequality. Assume this last inequality to be true. This would 
imply from Equation 3.7 and Equation 4.6 that 



N 



j Sg-(t) dt { J [ +q'(t) ] Sj-(t) dt } < — J 






T 



N. 



dt 



+ q*(t) 
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and rearranging this expression, we have 



[ j Sj“(t) dt ]- < J [ -2- + q-(t) ] Sj-(t) dt | 
T T - T 




dt 




q*(t) 



which is precisely a form of the Cauchy-Schwarz inequality applied to the inner 
product of two functions gj(t) and g-,(t). where in this case 



gi(0 = 

g.(t) = 



N 



J + q-(t)]'"‘ s.(t) 



s,(t) 



^0 



I + q"(t) 



The unsophisticated jammer case discussed in the previous chapter, where q(t) = 
K. in this case yields from Equations (4.3) - (4.5) 



NT 



4- K- 



E' = 



^ + K" 



P = 



Sn(t)s,(t) 



NT 



dt = p 



+ K‘ 



so that 



Pe,o = 



erfc [ 



E(1 -p) 



2(Nq + 2K*) 



] 



(4.7) 



which is identical to Equation 3.10. This is again not surprising since the assumed 
conditions result in simple AWGN interference for which the receiver of Figure 1.1 is 
indeed optimum, hence the result that Pe o “ ^e s‘ 



ANALYSIS OF RECEIVER PERFORMANCES WITH EXAMPLE 
INTERFERENCE MODELS 



The results of the two previous chapters are now used to evaluate receiver 
performance under various interference conditions for two important signaling 
schemes, namely Phase Reversal Keying (PRK) and Frequency Shift Keying (FSK). 
For PRK. the transmitted signals are 



Sj(t) = Am sin 2:tf^t - (-1)‘ A Vl - nr cos 2:tf^t 0 < t < T . i = 0.1 (5.1) 

where 0 ^ m < 1. The vector diagram of Figure 5.1 shows that the parameter m 
controls the phase angle between the signals Sp(t) and Sj(t). When m = 1. the two 
signals are undistinguishable. whereas for m = 0, the angle between the signals is ISO*^ 
and we have the familiar Phase Shift Keying (PSK) scheme. For FSK the transmitted 
signals are 



S;(t) = A sin 2nft 0 < t < T . i = 0.1. (5.2) 

We assume for simplicity that 2fT as well as 2fT (for i = 0,l) are integers. This 
means that 

T 

J Sj"(t) dt = .A'T ^ = E = Energy per bit i = 0,1 (5.3) 

0 

for both signaling schemes. Observe therefore that from Equation 1.3. p = 2m^ - 1 for 
PRK and p = 0 for PSK. 

Several interference conditions are now analyzed for the PRK and FSK signaling 
schemes introduced. For ease of comparison, all interferences are normalized so that 
they all have equal power (see Equation 2.4). Furthermore, we assume (as a worst 
case situation) that the interferer has timing information about the receiver (that is. the 
interference is synchronized to the receiver clock), and that it uses a waveform q(t) (see 
Equation 2.1 ) that repeats itself everv' T seconds. 
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COS27tf 



Figure 5.1 Vector Diagram of Signals Sjj(t) and Sj(t). 

A. LINEAR WAVE 
We define 



q(t) = 



{ 



a 

1 + b 



T 



0 



0 < t < cT 



cT < t < T 



0 <c < 1 



resulting in an interference generating waveform as depicted in Figure 
Equation 2.3 and Equation 2.4, we can obtain 

1 c 

Q^iO) = — J [ q(t) ]- dt = — (a-c^ + 3abc + 3b^) 

1 Q J 

for the average power of the nonstationar\’ interference. 



(5.4) 

2. From 
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Figure 5.2 Several Variations of Linear Wave Interference 
as a Function of (normalized) Time. 
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1. Phase Reversal Keying 

The linear wave interference of Equation 5.4 is now used to determine the 
receiver probability of error for PRK focusing first on the performance of the receiver 
of Figure 1.1. From Equation 5.1, we obtain 



Sj(t) = 2A yi - nr cos 2nf^t 0 ^ t ^ T. 

The following integral must first be evaluated in order to obtain namely 

cT ■’abt 

J Q(t)Sj-(t)dt = J (— + Ap- + b-) [ 4A-(1 - nr) cos' 2nf t ] dt 

T 0 

= 2A-T( 1 - nr) (a"c- + 3abc + 3b-)- 
3 

3(a"c-+ abc ) cos 47Tf cT - 3abc + 2{a"c- + 3abc+ 3b-){2irf cT)- 
2(a-c- + 3abc + 3b')( 2;rf cT)- 

3a"c- sin 4nf cT 3(a-c" + 2abc + b") sin 4;if cT 

r h 1 E 1 

4(aV+3abc+3b-)(27ri;cT)^ 2(aV+ 3abc + 3b') 2;rf cT 

{ 5 . 5 ) 

a. Performance of the Suboptimum Receiver 

From Equation 3.8 and Equation 5.5. we obtain 

JSR-SNR-c, 

C = ! (5.6) 

1 + JSR-SNR-Sj 

where 



JSR — (aV + 3abc+3b-).E 

E 3 



SNR = 



X ^ 

.Nq - 
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and 



*^1 = 



3{a“c"+ abc ) cos d/rf^cT - 3abc + 2(a"c* + 3abc+ 3b*)(27if cT)" 



3{a'c‘ + 3abc+ 3b")(2;ircT)' 



3a"C" sin 4k f cT 

r 



4{a-c- + 3abc + 3b- X 2Kf^cT/ 



3(a'c- + 2abc + b-) sin 4Kf cT 
H i C 1 

2(a-c- + 3abc + 3b‘) 27tf^cT 



This expression demonstrates that the performance of the conventional 
suboptimum receiver depends on the SNR as well as JSR and on the actual values of 
the parameters a. b. c. f^, and T. The receiver designer has limited control of JSR. 
Clearly, as JSR C 1 for constant SNR. and ^ -* 12 {see Equation 3.7). If 

the jammer is not present C = 0 and ^ = Pg- However, for constant JSR'SNR, the 
receiver designer can attempt to minimize C by proper choice of the product f cT. For 
fixed JSR. it is apparent that Ej 0 as f T ^ (for c constant). Thus C -+ 0 as f^T 

^ resulting in improved (suboptimum) receiver performance, namely P^ ^ P^. In 

practice, this means using as high a signaling frequency as possible, or choosing as long 
an observation time T as possible. 

Actual performance (in terms of probability of error) for the conventional 
receiver is obtained from Equation 3.7. Observe that for PRK. 



E(l-p) .A-T(l-nr) 



Nr 



Nr 



E 

— 5(1 - nr) = SNR(1 - nr) 

N'o 2 



(5.7) 



and from Equation 5.6 



1 - C 



1 

1 + JSR-SNR-e, 



so that 



P(, g = erfc ( v^SNR{ 1 - m')a 2 ] 



(5.S) 
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probability of error 

4 dT • 4 “ 4 ^ 4 • 



where 



a = I (I + JSR-SNR-Cj). 

A plot of Pg j as a function of SNR and JSR is shown in Figure 5.3 with m set to zero 
(resulting in PSK modulation), and parameters b set to zero and c set to one. 
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Fiaure 5.3 Performance of the Suboptimum Receiver for PSK .VIodulation 

with Linear Wave Interference. 



b. Performance of the Optimum Receiver 

The previous result for P_ . can be compared to the performance of the 
optimum receiver. From Equation 4.5 and Equation 4.6, we have 
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E'fl-p) A-Tfl - nr) 

-ij 



1 + cos 4niyi'x 



Nr 



1 + (a"X“ + 2abx + b") — ^ 



dx 



1 

+ J ( 1 + cos dJlf^Tx) dx ] 
c 



(5.9) 



where the change of variables x = t T has been used in Equation 5.9. If we let b = 0 
and c = 1 as a special case, then q(t) = at;T, 0 ^ t ^ T, Q^IO) = a' 3, and JSR = 
(a" 3) E so that 

E'(l - p') A-T(l - nr) } 1 + cos 4rtfTx 

= ^ J dx. 

- ^'o 0 1 + 3JSR-SNR x- 



Let v^jJSR'SNR x = y. then dx = dv4/3JSR*SNR, so that 

V3JSR-SXR 

E'(l - p') A"T(1 - nr) 1 , 1 + cos (dttf Tv y3JSR-SNR) 

= : — j '' ' dy 

2 N y3JSR'SNR 0 \ + y~ 

which can be expressed in the form 



y3JSR-SNR 

E'(l-P') SNR(l-nr) , , cos (4jif Tv 

^ f tan'* V3JSR-SNR + J ^ 



/ 



3JSR-SNR) 



V3JSR-SNR 



1 + v- 



dy] 



(5.10) 



where the remaining integral in Equation 5.10 can only be evaluated numerically or an 
approximation can be obtained, if it is assumed that y3JSR*SNR » 1. (The derivation 
of such an approximation is presented in Appendix A). In the case of JSR = 0. then 
q(t) = 0. so that directly from Equation 4.6 we have 
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H'( 1 - p ) 



Ed - p) 



(5.11) 



2N 



J Sj-(t) dt = SNR(I - nr) = 



0 T 



NV 



Thus we obtain 



erfc [ 7 SNR(1 - mf) 2 ] JSR = 0 

erfc [ V'SNR(1 - nr)P 2 ] JSR s 0 



where 



P = 



73JSR-SNR 



[tan'V3JSR-SNR + J 



73JSR-SXR 

cos (4JtfTy 73JSR-SNR) 



1 + V- 



dv 



* 7t 73JSR-SNR . 



A plot of Pg Q as a function of SNR and JSR is shown in Figure 5.4 with m set to 
zero, and parameters b set to zero and c set to one. 

Comparison of Equation 5.8 and Equation 5.12 is dilTicult at best because 
the parameter p can only be approximated (see .Appendix A). We observe that with b 
= 0 and c = 1, 



a = 



1 + JSR-SNR- 



3 + 2(2;rfJ)‘ 

2(2;ti;T)- 



If JSR is very large (with SNR fixed) and iyr is allowed to become 
unbounded (in order to the suboptimum receiver to overcome some of the jamming 
present), then 



a = 1 (1 + JSR-SNR) 
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while 



P 



K 

v'3JSR-SNR 



is a close approximation, and clearly P > a. 




Figure 5.4 Performance of the Optimum Receiver for PSK Modulation 
with Linear wave Interference. 



For example, when JSR’SNR takes on a value of 15dB, P a ~ 10, 
demonstrating the fact that P^ ^ < P^ ^ (as expected). One must however note that as 
JSR increase without bound, for fixed SNR. both a and P tend to zero resulting in 
both receivers displaying an increasingly worse performance, namely Pg ^ = Pg s 
0.5. 
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2. Frequency Shift Keying 
From Equation 5.2, we have 



Sj(t) = 2 A cos 7t(fj + Fq)! sin 7t(fj - f^it 



0 < t < T 



and we evaluate 



cT 

I’ Q(t)Sj-(t)dt = f q-(t) [ 2A cos trffj + rp)t sin 7t(fj - fg)t ]- dt 
' T ■() 



.A'Tc 

= — ; — (a‘c‘ + 3abc+ 3b“)* 
3 



3{a"c + ab) 2cos Ttfaj-ttplc cos2;ta|C cos27taQC 2cos + ttglc 
; — ^ ^ 1 ^ ^ ^ 



7i^c(a^c^ + 3abc-h 3b^)^ 


+ 

(«i-V 


v+' 

(2a,)- 




■ i U j 

(«1 + «q)- 


3(a-c* + 2abc + b") 

r 


2sin 7t(a,-a,,)c ^ 


sin27ta,c ^ sin27tapC 


2sin 7t(a, -t-a^ic 


2;ic(a'c^-h 3abc+ 3b^/ 


«r«o 


2a, 


2aQ 


«l+«0 


3a^ 

r 


2sin 7t(a,-a,,)c ^ 


sin2;ta,c 

1 J_ 


sin2;taQC 


2sin 7t(a, +a^)c 


7t- c(a"c- + 3abc + 3b“)^ 




(2a, 


(2a„)-^ ■ 


(aj + ao)' 


3ab ^ 


2 


1 


1 


■> 


, S S . ^ . . .. S [ 


. s ^ 




- s 


s I 



7T-c(a'c- + 3abc+ 3b*)^ (a^-aQ)- (2a,)- (2ag)- (aj + ttp)- 



(5.13) 



where a, = 2f,T. i = 0,1. 

a. Performance of the Suboptimum Receiver 

From Equation 3.8 and Equation 5.13 we obtain 

JSR-SNR- 

C = ! — 

1 H- JSR-SNR-c,' 



( 5 . 14 ) 



where as before 



JSR = [-^(a-c- + 3abc+3b-)] E 



SNR = 



^0 - 



and 



c,' = (I 



3(a‘c + ab) 2cos n(a|-ap)c cos27rajC cos2n:apC 2cos ;r(aj + af,)c 
K-c(a-c' + 3abc + 3b') (aj-ag)' (2aj)‘ ^ 



3(a'c' + 2abc + b') 2sin iT(a,-aQ)c sin2Jra|C sin2n:a^c 
TT~~r~. , ^ ' 



(«i + «o)‘ 

2sin K(a| -t-a^lc 



2Kc(a'c* + 3abc+ 3b') 



0 



2a 



2a 



0 



«l+«0 



3a- 






7T'^c(a'c' + 3abc + 3b') (aj-Og)'^ 

3ab 2 



sin 3T(aj-ag)c^sin2Ka|C ^sin27ragC 2sin irfaj +«g)c 



(2a, (2a,)^ 



(ai+ao) 



. 4 -- 



• 4 -- 



7t‘'c(a“c" + 3abc+ 3b“) (a^-ag)"' 



/ 1 



In order to maximize C, Cj' must be minimized, with the previously 
introduced restriction on and that they be integers. For b = 0, and c = K 



2cos ;r(a,-an)c cos2;ia,c cos2;ranC 2cos ;r(ai 
c,' = 1 - 3[ -— 1- 1 -Sl I - + ^ ^ il— ] (5.15) 

' K-(a,-ag)' (2;raj)- (27rag)' 7r-(aj+ag)' 

so that minimization of Cj' can be accomplished by maximizing the term in brackets in 
Equation 5.15. It appears that optimum choices for these parameters are values that 
result in a^ - Ug being small. Table 1 shows that when ttj = 3, ttg = 1. c = 1 and b = 
0. then C|' = 0.801579349, the minimum value achievable. Thus the optimum value of 
C at li.xed J SR- SNR is 



0.8016 JSR-SNR 

C = (5.16) 

1 + 0.8016 JSR-SNR 

which can be achieved with finite values of f^ and T. unlike the PRK scheme, where 
infinite frequency-time products are required for optimum results. 
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TABLE 1 



COMPUTED VALUES OF c,' FOR GIVEN VALUES OF a, AND a 



= 



«i -«o> 



■•0 



«1 








2 


1.4454 


3 


1.5562 


4 


1.5823 


5 


1.5926 


6 


1.5978 


7 


1.6007 


8 


1.6025 


9 


1.6037 


10 


1.6045 





"l' 




(«, = 3) 


4 


0.9625 


5 


1.0331 


6 


1.0495 


7 


1.0562 


8 


1.0597 


9 


1.0618 


10 


1.0631 


11 


1.0640 


12 


1.0647 



«1 






II 




0.8016 


4 


0.8412 


5 


0.8460 


6 


0.8472 


7 


0.8476 


8 


0.8478 


9 


0.8479 


10 


0.8479 


11 


0.8480 



«1 


II 


5 


0.8999 


6 


0.9054 


7 


0.9581 


8 


0.9603 


9 


0.9611 


10 


0.9615 


11 


0.9617 


12 


0.9618 


13 


0.9619 



«d 


^l' 

(ao=l) 




«d 


Q 

^ n 

11 


5 


0.9338 




10 


0.9215 


6 


0.9151 




11 


0.9249 


7 


0.9277 




12 


0.9224 


8 


0.9197 




13 


0.9245 


9 


0.9257 
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Actual performance (in terms of probability of error) for the conventional 
receiver is obtained from Equation 3.7. Observe that for FSK 



E(1 - pi 






= SNR 2 



(5.17) 



and 



1 

1 - C = 

1 + JSR-SNR-C,' 

so that 

Pg ,, = erfc ( v’SNR a 4 ) (5. IS) 

where 

a = 1(1 + JSR-SNR-Cj). 

■A plot of Pg j as a function of SNR and JSR is shown in Figure 5.5 with parameters b 
set to zero, and c set to one. 

b. Petfortnance of the Optimum Receiver 

Again, the above result for the performance of the suboptimum receiver 
can be compared to the performance of the optimum receiver using Equation 4.5 and 
Equation 4.6. For the FSK modulation case being considered and linear wave 
interference. 

E'(l 



T 

+ J [ 2.A cos 7T(fj + fpit sin 7t(fj - IqH ]" dt }. 
cT 



p') 1 [2.A cos 7r(fj + fp)t sin Ttifj-fplt]' 



2N4 






1 + (-^t- + 2abt + b-);^ 
T" 2 



dt 
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Figure 5.5 Performance of the Suboptimum Receiver for FSK Modulation 

with Linear Wave Interference. 

This expression can be put in compact form with b set to zero and c set to one as a 
special case. Some algebraic manipulations yield 



E'(l - p')_ ^SNR ^ 



V3JSR-SNR 



4 j. [cos (t 1«[ + aQ)x/y3JSR’SNR) sin (-t(«j-«q)x;V3JSR'SNR) ]“ ^ ^ 

J ^ T~ — i ^ ■ 

1 + X" 



V3JSR-SNR Q 



(5.19) 
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This equation can only be evaluated numerically for set values of the 
parameters and and changing values ot'JSR’SNR. When = 3 and Og = 1 
(which is the optimum choice for minimizing C and thus minimizing P^, ^ ). Equation 
5.19 becomes 



y3JSR-SNR 

E'(l-p') 2SNR , [cos (2rtx y3JSR-SNR) sin (rtx 73JSR-SNR) ]“ 

J ^—5 dx]. 



v'3JSR*SNR 



1 + X- 



(5.20) 



[f JSR = 0. then q(t) = 0. therefore E'(l -p');2 = E(1 - p)'Ng. Thus 



E'(l - p ) 
•) 



E 2 






SNR 2 



(5.21) 



so that 



erfc ( y SNR 4 ) JSR = 0 

^ erfc ( y^sR p 4 ) JSR s 0 



where 



P = 



V 3JSR-SNR 



y3JSR-SNR 

, [cos (2rtv v^JSR-SNR) sin (Ttv V3JSR-SNR) ]- 
J ^ ^ dy. 



1 -h V" 



A plot of 0 ^ function of SNR and JSR is shown in Figure 5.6 with parameters b 

set to zero, and c set to one. Results of numerical evaluation of P are presented in 
Table 2 for b = 0 and c = 1. and compared with the values a. 
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Fisure 5.6 Performance of the Optimum Receiver for FSK Modulation 
with Linear Wave Interference. 



3. Comparison 

Performance comparison of the receivers for PRK and FSK with linear wave 
interference, from Equation 5.6 and Equation 5.16 shows that the parameter C for 
FSK does not grow as rapidly to 1 as a function of JSR'SNR. as does C for PRK. 
Thus, the frequency diversity of FSK provides a small jamming margin over PRK. 

We note that dilTerences between the parameters P and a are much smaller for 
FSK than dilTerences between the corresponding factors P and a in the PRK scheme. 
Hence suboptimality is somewhat reduced in the FSK scheme. 
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TABLE 2 

COMPARISON OF VALUES OF a AND p FOR GIVEN VALUES OF JSR-SNR 



JSR-SNR 




2 


10 


50 


100 


200 


1.000 


a 




0.3841 


0.1109 


0.0243 


0.0123 


0.0062 


0.0012 


p 




0.4214 


0.1435 


0.0387 


0.0212 


0.0114 


0.0026 


Pa 




1.0970 


1.2938 


1.5898 


1.7205 


1.8390 


2.0867 



B. SINUSOIDAL WAVE 

We define 



q(t) = 




yabc 



+ b- 



0 < t < cT 

0 < c < 1 

cT < t <T 



(5.23) 



as another interference generating function which can be depicted as shown in Figure 
5.7. Then from Equation 2.3 and Equation 5.23 we obtain 

1 "T 
- — Jo f 

= — ; — ( a“c" + 3abc + 3b" + a"c" sine 4f^cT + 2acV6(abc + b“) sine 2f^cT ] 



for the average power of the nonstationar\' interference. Observe that whenever 2f,cT is 
an integer. Qg(0) above is identical to the average power of the linear wave 
interference. 
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Figure 5.7 Several Variations of Sinusoidal Wave Interference 
as a Function of (normalized) Time. 
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1. Phase Reversal Keying 

For PRK. from Equation 5.1 we obtain again 



Sj(t) = 2A yi - m" cos27tf t 0 ^ t ^ T. 



Prior to obtaining receiver performance we evaluate the integral 
cT 

J Q(t)Sj‘(t)dt = J q"(t) 4A“(l-m") cos“ 2rtf t dt 
0 



= 2A“T( 1-m") — [ a“C' + 3abc+ 3b‘ + a‘c“sinc 4f cT + 2acv^6(abc+ b~) sine 2f cT ]* 

3 c t 

0.5a"C"[ 2sinc 4f cT + sinc4(f-f)cT + sine 4(f + f )cT ] 

f I _j_ T c r c r * 

a^c" + 3abc + 3b" + a"c^ sine 4f^cT + 2acV6{abc+ b") sine 2f^cT 



acy6(abc+b") [ sine 2(f^-2gcT + sine 2(f^ + 2gcT ] 
a"c" + 3abc + 3b" + a"c" sine 4f^eT + 2acy6(abc + b^) sine 2f^cT 

3(abc + b") sine 4f cT 

+ _ }• 

a"c" +3abc +3b^ +a"c" sine 4f^cT + 2acv o(abc+ b") sine 2f^eT 



where 



(5.24) 



sin Jtx 

sine (x) = 

;tx 



a. Performance of the Suboptimum Receiver 

From Equation 3.S and Equation 5.23, we obtain 



JSR-SNR-c 

C = ^ (5.25) 

1 + JSR-SNR-e 

S 
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JSR = 



— [ a*c“+ 3abc+ 3b* + a*c* sine 4f^.cT + 2acV6(abc + b") sine 2f^,eT ] 



E 



SNR = 



No2 



and 

0.5a*e*[ 2sine 4f eT + sine 4(f -f)cT + sine 4(f +f)eT ] 

+ 1 £ l.sL_£ r 

a*e* + 3abe + 3b^ +a*e^ sine 4f^eT + 2ae^/6(abe + b*) sine 2f^eT 

aev/6(abe + b“) [ sine 2(f^,-fj.)eT + sine 2(f^+ f )eT ] 
a*e* + 3abc + 3b* + a*e* sine 4f^eT + 2aev^6(abe+ b") sine 2f^eT 

3(abe + b*) sine 4f eT 

+ ^ , }• 

a*e* + 3abe + 3b* +a*e* sine 4f^eT + 2aev'6(abe + b*) sine 2f^eT 

(5.26) 

For the speeial ease in whieh b = 0 and e = 1, then = 1 for ?= f and 
C = 1.5 for f = f . thus we have 

s c r 

/ JSR-SNR (1 + JSR'SNR) if 

I 

( 3JSR-SNR (2 + 3JSR-SNR) iff^ = f^. 

It is apparent that C is largest when f^ = for fi.xed JSR*SNR. This means (as eould 
be expeeted) that the jammer is most damaging to the reeeiver when it operates at the 
signal frequeney, and furthermore there is nothing the suboptimum reeeiver ean do (in 
terms of ehoosing longer integration times for instanee) in order to reduee jamming 
elTeets, other than to switeh to a new frequeney of operation. 



Actual receiver performance is obtained from Equation 3.7. Observe that 



for PRK 




(5.27) 



and 



1 - C = 1(1 + JSR-SNR-c ) 

^ S 



so that 



P 



erfc [ ySNR(l - nr )a 2 ] 



(5.28) 



where 



a = 1 (1 + JSR-SNR-c^). 

.A. plot of Pg j as a function of SNR and JSR is shown in Figure 5.8 with m set to zero, 
f^ = f and the parameters b set to zero, and c set to one. 
b. Performance of the Optimum Receiver 

Again, we compare previous result with the performance of the optimum 
receiver using Equation 4.5 and Equation 4.6. We obtain 



E'(l - p ) 



1 . 4A"( 1-nr) COS" 2ttf t 

— 1 ^ — dt 




cT 



4.A"( 1-nt") cos" 27tf t 

^ r 





dt 



T 

+ j [ 4A"(1 - nr) COS" 2Jtf t ] dt }. 
cT 
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Figure 5.8 Performance of the Suboptimum Receiver for PSK Modulation 
with Sinusoidal Wave Interference. 



If we assume that b = 0 and c = 1 as a special case, then q(t) = cos 2a:f^t, 

QJO) = a^3. and JSR = (a-;3), E so that 



E'(l -p') 1 



— = — f 



4A"( 1-m") COS" 2;if t 

r 



0 0 



N 



1 + (1 +COS 4;tf t), — ^ 



0 



A-{l-m-)n 

2(a-3) 




0 



1 + cos 4;tf t 

dt 

1 + n cos 4Jtf t 

* C 



(5.29) 
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where 



a- 3 _ JSR-SNR 

a" 3 + Nq 2 ~ 1 + JSR-SNR 



0 < n ^ 1- 



In order to evaluate the integral in closed form, we use the tabulated integral 
[Ref. 3; pp. 107-122], 



1 cos q dx 

2rc Q I + ri cos x 






1 + P- 

1 - p2 



P = 



1 + n/i - 



0 < n < 1. (5.30) 



The integral in Equation 5.29 can be put in the form 



J 

0 



T 



1 + cos 4;tf^t 
1 + n cos 4;tf t 

* C 



4;tf T 
dt=J ^ 



0 



1 + cos qx 
1 + r\ cos X 





In order to obtain a closed form expression for this integral, we assume 
that q is an integer (The most interesting case is q = 1, so the restriction is not too 
unreasonable). With this assumption, the integrand is periodic of period 2n, and since 
4;tf T = 2n'2f T and 2f T is assumed to be an integer, sav 2f T = n, then 

c C C ^ c 



,T 1 + cos 4;tf t 
Q 1 + n cos 4nf^t 



dt 



1 

4;tf 

C 




1 + cos qx 

^ dx 

1 + r| cos X 



T n 

iiVr'77 



2n 



1 + cos qx 
1 + x\ cos X 



dx = T[ 1 +(-p)^ 



1 - p- 



and finally 



E'(l -p') 



= SNR(1 - m2)(l 1 + (-p)fi 



1 + P^ 

1 - p‘ 



= SNR(1 -m“)( 1 + (- 



-n .n ■ n)(i + - n' ) 



r)"! 



1 + yTtp" 1 - n* +Vi - 



= SNR(1 - m-) V 



(5.31) 
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where 



1 + J 1 + 2JSR-SNR 

V = — 

(1 + 2JSR-SXR) + (1 + JSR-SNR) 71 + 2JSR-SNR 



If JSR = 0 then q(t) = 0, so that directly from Equation 4.6 we have 



E'(l-p') A-T , 

= ( 1 - irr). 

■) \ 

.Nq 



(5.32) 



Thus we obtain probability of error for the optimum receiver, namely 



erfc [ 7 SNR(1 - m^ 2 ] 
erfc [ 7 sNR( 1 - m^ 2 ] 



JSR = 0 



JSR * 0 



(5.33) 



where v is defined following Equation 5.31 A plot of ^ as a function of SNR and 
JSR is shown in Figure 5.9 with m set to zero, and the parameters b set to zero, with c 
set to one. 

If we compare the parameters v and a when f^. = f^ (i.e.. q = 1). where b 
= 0. and c = 1, we find that 

2 1 + 7l + 2JSR-SNR 

a = < — 

2 + 3JSR-SNR (1 + 2JSR-SNR) + (1 + JSR-SNR) 7l + 2JSR-SNR 



(5.34) 

Observe that this inequality becomes stronger as JSR increases with 
constant SNR. For instance at a JSR-SNR of 15dB. the two sides differ by a factor of 
1.34. and at a JSR-SNR of 25dB. the factor becomes 1.44. Nevertheless, both sides of 
the inequality tend to zero as JSR co. implying both ^ and Pg q under 

that conditions. The optimum receiver has a definite "advantage" over suboptimum 
receiver for high JSR . however this "advantage" has a limit as v a 1.5 as JSR -» ^ 
in Equation 5.34. 
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2. Frequency Shift Keying 

The above result for the performance of the suboptimum receiver can be 
compared to the performance of the optimum receiver using Equation 4.5 and 
Equation 4.6. For the FSK, from Equation 5.2 we have 

Sj(t) = 2A cos 7t(fj + fg)t sin 7T(fj - fg)t 0 < t ^ T 
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and first evaluate 



cT 

J Q(t)Sj"(t) dt = J q“{t)[ 2A cos 7t(fj + fg)t sin 7t(fj - fg)t ]" dt 

j 0 

= A"T( — [ a'c" + 3abc + 3b" + a"c" sine 4f' cT + 2acv^6(abc+ b") sine 2f cT] }• 

2, c c 

fj 4 - + 3b")(2sinc 2f^cT - 2sinc 2fjcT + sine dfjCT + sine 4fgcT) 

2(a*c" + 3abc + 3b" +a"c" sine 4f^cT + 2acy6(abc + b") sine 2f^cT) 



+ a-c-[- 



sine 2(2f^-f^)cT + sine 2(2f, + f^)cT - sine 2(2fArg)cT - sine 2(2f^,4-fj)cT 
2(a"c" + 3abc +3b" +a\’" sine 4f^cT + 2acv^6(abc+ b") sine 2f^eT) 



sine 4(f^-f|)cT + sine 4(f^+ fj)cT + sine 4(r,-fg)cT + sine 4(f^ + fg)cT 
4(a"c" +3abc +3b" +a"c" sine 4f^cT + 2acV6(tibc+ b") sine 2f^cT) 



1 — sine 2(f -f )cT + sinc 2(f, + r,)cT-sinc 2(f,-f. )cT-sinc 2(f.4-f )cT 

+ 2ac76(abc4-b")[ ^ ^ ^ ^ 

2(a"c" + 3abc + 3b" + a"C"sinc4f^cT + 2aeV^6(abe + b")sinc2f^cT) 



sine 2(f^.-2fj)cT + sine 2{f^+2f|)cT +sinc 2(f^-2fg)cT + sine 2(r^+2fg)cT 
4(a"C" + 3abc +3b" +a"c" sine 4f^cT + 2ae v^6( abe + b") sine 2f^,cT) 



(5.35) 



where = fj + fg and fg = - fg- 



a. Performance of the Suboptimum Receiver 
From Equation 3.S we obtain 



JSR-SXR-c ' 

C = ^ (5.36) 

1 + JSR-SNR-e ' 
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where 



JSR = 



— [ a‘c" + 3abc+ 3b' + a'c' sine 4f cT + 2acv'6(abc + b') sine 2f eT 
3 c c 



SNR = E (N- 2) 



and 



£s = 1 + 



(a'e" + 3abe+ 3b')(2sine 2f^eT - 2sine 2fjeT + sine 4f|eT + sine 4fpeT) 
2(a'e" + 3abe + 3b" +a'e" sine 4f eT + 2ac'v/6(abe + b") sine 2f eT) 



-) ^ sine 2(2f -f )eT + sine 2(2f +f )cT - sine 2(2r-f,)eT - sine 2(2f + f )eT 

I C ^ C ^ C q ^ C a' 

^ ^ ^ L ^ ^ 

2(a"e" + 3abe + 3b^ +a'e^ sine 4f^eT + 2aeV6(abe + b") sine 2f^eT) 



+ 



sine 4(f.-f|)eT + sine 4(f, + fj)eT - sine 4(f,-fQ)eT + sine 4{f^ + f(^)eT 
4(a'e" + 3abe + 3b' +a'e" sine 4f eT + 2ae>/6(abe + b') sine 2f,eT) 



j =r sine 2(f -f )eT + sine 2(f +f )eT-sine 2(f -f',)eT-sine 2(f +f.)eT 

+ 2aev6{abe+b-)[ ^ — ± — 

2{a"e^ + 3abe + 3b" + a^e^sine4f,eT + 2aey6(abe+ b')sine2f^eT) 



^ sine 2(f,-2f| )eT + sine 2(f,+ 2f|)eT + sine 2(f,-2fQ)eT + sine 2(f^+ 2f^)eT 
4(a'e' + 3abe +3b' +a'e' sine 4f^eT + 2aev'^6(abe-^ b') sine 2f^,eT) 

(5.37) 



Note that C in Equation 5.36 ean be ntinimized by setting c^' to its smallest 
possible value. Sinee with b = 0 and e = 1 £^' equals 1.5 For 2, 0.5 for = 

fj 2. and 0.75 for f^ = fj or f^ = f^. it is apparent that the januner is least damaging 
when f^, is set to (f^ - Fq) 2, and most damaging when f^ = (Fq + Fj ) 2. the midpoint of 
the two operating Frequeneies. From Equation 3.7, we obtain aetual reeeiver 
perFormanee. Observe that For FSK 
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PROBABILITY OF ERROR 

4 nT® 4 4 a 



E( 1 - p) 



X 



0 



SXR 



so chat 



{5.3S) 



= erfc ( ysXR a, 4 ) 



(5.39) 



where a = 1/(1 + JSR*SXR*Cj ). A plot of ^ as a function of SNR and JSR is 
shown in Figure 5.10 with = SCq, and the parameters b set to zero, and c set to one. 
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h. Peyfonnance of the Optimum Receiver 

The result of Hquation 5.39 can once again be compared to the 
performance of the corresponding optimum receiver using Equation 4.5 and Equation 
4. 6. We thus obtain 



E'(l-p) 1 



2N 



Hi 



cT 



2A cos n(fj + fp)t sin 7t{f|-fj,)t ]“ 



0 0 



1 + ( /— ac cos 2rtf^t + yabc+ b' ) 






dt 



T 

+ J [ 2.A cos Jt(f^ + fp)t sin Tt{fj - fp)t ]" dt } 
cT 



If set c to one and b to zero as a special case, then 



E'(l - p’) 1 [ 2A cos 7t(f, + fp)t sin K(f, -f,)t ]* 



2N 



j ‘ dt 



0 0 



I + (— + —cos 47tf t) 

4 ^ C ^ 



A~\\ c [ 1 + COS (f^ X 2f^) ] sin* (f^ x 4f,) 

^ ■fUTTITTr 3 ^ 



4ttf^(a- 3) 



1 + 11 cos X 



dx 



(5.40) 



where as before rj = JSR’SNR (1 + JSR'SXR). and 2f^T is assumed to be an integer, 
say 2f^T = n. Tlie integral of Equation 5.40 must be evaluated numerically or it can be 
approximated under certain circumstances. (The derivation of such an approximation is 
presented in Appendix B }. If JSR = 0, q(t) = 0, then 

E'(l-p') SNR 



so that the actual receiver performance is obtained from 



erfc ( y SNR 4 ) JSR = 0 




erfc ( ySNR p 4 ) JSR » 0 
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PROBABILITY OF ERROR 

i nP® « W”® « 4 rV® 



where 



[ 1 + COS (f^ y 2f.) ] sin- (f^ y;4f ) 
1 + t] cos y 



47if T 



p = f 

2;ifJ(l+JSR-SNR) 5 



dy 



2(1 + JSR-SNR) + 27i + 2JSR-SNR 
2 + 5JSR-SNR + 2 (JSR-SNR)^ + (2 + 3JSR-SNR)Vl + 2 JSR-SNR 



A plot of Pg Q as a function of SNR and JSR is shown in Figure 5.11 with parameters 
b set to zero, and c set to one. 




Figure 5.11 Performance of the Optimum Receiver for FSK Modulation 
with Sinusoidal wave Interference. 
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It is not dillicult to show that for f = f 2 and c = I is alwavs exceeded 

C S 

by p. However, for increasing JSR*SXR. the ratio P a 1.5. 

3. Comparison 

.Again, it becomes apparent that interference effects for FSK are not as severe 
as for PRK as the differences in performance between the optimum and suboptimum 
receiver for the former modulation scheme are not as large as those the latter 
modulation scheme. 



C. PULSED WAVE 

We define 



q(t) = 




4 - 2abc + b" 



b 

0 



n <t<(C+!/ 2 )T 
p ' p 

{e+‘y)Tp< t < {£+i)Tp 



0 < t < cT 



cT < t 1 

(5.43) 



as yet another form of an interference generating function where f = 0.1, p-1. pTp 

= cT. and 0 < c < 1. 

In this example, q(t) is being pulsed many times between two levels in the 
interval 0 ^ t ^ cT as shown in Figure 5.12. From Equation 2.3 



1 "T 

Q.,(0) = — ]■ q‘(t) dt 
* 0 




p-1 

f = 0 



2a~c~ 



+ 2abc +b-) 



(t+ '^2)Tp 



-I- 




P -1 
V 

f = 0 






(£+ 1)T, 



{£+ H)T 



dt 



-(a'C" + 3abc + 3b*) 



for the average power of the nonstationarv' interference. 
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Figure 5.12 Several Variations of Pulsed Wave Interference 
as a Function of (normalized) Time. 
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1. Phase Reversal Keying 

The Following integral must first be evaluated in order to obtain Pg namely 
cT 

]■ Q(t) Sj-{t) dt = J q-(t) [ 4A*(1 - nr) cos* 27tFt ] dt 
T 0 

= 2A"T{ 1-nr) — (a*c‘ + 3abc+ 3b*)‘ 



{1+ [sin 47tf cT + 2sin- 2;tf cT tan'^ (kF cT p)]j (5.44) 

4;rf;cT " 



(See Appendix C For the pertinent derivations). 

a. Performance of the Supoptimuni Receiver 

The pulsed wave interference of Equation 5.43 is now used to determine the 
receiver probability of error For PRK Focusing First on the performance of the receiver 
of Figure l.l. From Equation 3.8 and Equation 5.44, we obtain 



JSR-SNR-c 

C = Ei (5.45) 

1 + JSR-SNR-c 

p 



where 



JSR 



— (a'c* + 3abc + 3b') 
3 



] E 



SNR = 



No 2 



c = ( 1 ^ [ sin 4;tF cT + 2sin* 2;tF cT tan (ttf cT p) ] } . 

r 



56 



Actual receiver performance (in terms of probability of error) for the conventional 
receiver is obtained from Equation 3.7. Observe that for PRK. 



E( 1 - p) A^T 



NV 



N. 



(1 - nr) = SXR(1 - nr) 



(5.46) 



and 



1 - C = 



1 

1 + JSR-SNR-C 

P 



so that 



Pg 5 = erfc [ ySNR( 1 - m')a. 2 ] 



(5.47) 



where 



a = 1(1 + JSR-SNR-Cp). 



A plot of Pg g as a function of SXR and JSR is shown in Figure 5.13 with m set to 
zero, and the parameters b set to zero, and c set to one. 

b. Performance of the Optimum Receiver 

The previous result for P^ ^ can be compared to the performance of the 
corresponding optimum receiver. From Equation 4.5 and Equation 4.6, we obtain 



E'(l - p') 4A-(1 - nr) cos- 2Jtf t 

- - ^ ^dt 



2X. 



1 + q-(t) 



Xr 



j COS* 27tf t dt 
cT 



When b is set to zero and c set to one as a special case 



. then q(t) = V2a'. 3, so that 



57 



E'n - p'l 

*) 



A-T(l -nr) 



1 



a 



j 





) 



= SNR(1 -m-)( 



1 + JSR-SNR 
1 + 2JSR-SNR 




SND IN DB 



Figure 5.13 Performance of the Suboptimum Receiver of PSK Modulation 
with the Pulsed Wave Interference. 



In the case of JSR = 0. q(t) = 0, so that directly from Equation 4.6 we have 



E'(l - P') 1 , , E(1 - p) 

— = f s.-(t)dt = — =SXR(l-m-) 

-.Nq J .Nq 



(5.4S) 
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so that 



P 



e.o 



erfc [ y SNR(1 - m') 2 ] 
erfc [ ySNR(l - nr)P 2 ] 



JSR = 0 
JSR 0 



(5.49) 



where 



1 + JSR-SNR 

P = • 

1 + 2JSR-SNR 

A plot of Pg Q as a function of SNR and JSR is shown in Figure 5.14 with m set to 
zero, and parameters b set to zero, with c set to one. 

From previous results we observe that P^ ^ ^ s 

1 1 + JSR-SNR 

a = < =6 

1 + JSR-SNR 1 +2JSR-SNR 

where c = 1 has been used in the above e.xpression for a. Note that the above 
inequality becomes stronger as JSR-SNR increases. 

2. Frequency Shift Keying 

In order to obtain P^ ^ for FSK modulation, the following integral must 
evaluated, namely 



J Q(t) Sj-(t) dt = X 

T £=0 ^ 



(C+ ‘4)Tp 

+ 2abc+ 2b") j [4A“ cos" 7T(fj + fg)t sin" K(fj-fg)t] dt 
£T 

P 




P -1 

(aV+3abc+3by Y 
e=o 



(£+ >'.)Tp 



J ( 1 + 

£Tp 



2JtfTt 2;rfTt 

cos ^ cos ^ — 

T T 



1 4nf,Tt 1 4JtFTt 

— cos * cos ) dt. 

2 T 2 T 
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PROBABILITY OF ERROR 

4 /V® 4 4 < /N"** 




Figure 5.14 Performance of the Optimum Receiver for PSK Modulation 
with Pulsed Wave Interference. 



.A.11 integrals of cosine terms in the previous expression are of the form 



JTkt sin Jtkt.T (6+ W)T 

(T T tk,T Vt 

P P 



sin rtkc(E+ p - sin )rkc£. p 

itlr 



where k = 2fT. Therefore The sum 
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p -1 

^ [ sin 7rkc(C+ ‘ p - sin (TtkcC p) 

f = 0 



(5.50) 



, rtkc Tike 

= sin" tan o 2sin ttke 

2 4p 



(See Appendix D for the pertinent derivations). Therefore 



A"Tc 

I Q(t) Sj"(t) dt = (a*c* + 3abc +3b")* 

T ^ 



1 , 7r(f, + f)cT 

[ sin- 7t(f, + f„)cT tan — ^ + 2sin 2rt(f, + f„)cT ] 

jr(f, + fQ)cT ‘0 2p W 0' J 



1 , ;r(f, - f,,)cT 

— - — — — [ siir 7t(fj - fQ)cT tan — ^ — + 2sin 2n:(q - 1 q)cT ] 

7t(fj - fQ)cT ' ^ 2p ‘ ^ 



1 Tif.cT 

( sin" 2;tf, cT tan o 2sin dirf.cT ) 

2jrfjCT '■ p '■ 



2^f'oCT 



( sin" 27rfQcT tan 



nf„cT 

— 2 h 2sin 47 t(qcT )) . (5.5 f) 



a. Performance of the Suboptimum Receiver 

From Equation 3.S and Equation 5.51, we obtain 



C = 



JSR-SNR-s ' 

n_ 



1 + JSR-SXR-e_ 



(5.5 , 



where 



JSR = [ — (aV + 3abc + 3b-) ] E 
3 



SNR = 



E 



> 0 ^ 
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and 






1 + 



7t(f, + fo)cT 



[ sin" 7T(f, + fn)cT tan — * — 1- 2sin 2;t{f, + fn)cT ] 

Jt(fj + fo)cT ‘ ° 2p Mo') 

1 7t{f, - L)cT 

[ sin- 7t(fj - fy)cT tan — ^ + 2sin 2K{fj - f^KT ] 



R(fj- fpicT 
1 

27tfjcT 



27tfQcT 



-P 



Jtf.cT 

( sin" 2JtfjcT tan + 2sin 4;tr,cT ) 



M' 



, ttLcT 

( sin" 2;tfQcT tan 1- 2sin d^f^cT ). 



(5.53) 



Actual receiver performance (in terms of probability of error) is obtained from 
Equation 3.7. Observe for FSK modulation 

E(l-p) SNR 



and 



1 - C = 



I 

1 + JSR-SNR-c ' 
P 



so that 



Pg = erfc ( VSNR a 4 ) 
where 



(5.55) 



a = 1(1+ JSR-SNR-c '). 

p ' 
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A plot of Pj, j as a function of SNR and JSR is shown in Figure 5.15 with parameters 
b set to zero, and c set to one. 




Figure 5.15 Performance of the Suboptimum Receiver for FSK Modulation 

with Pulsed Wave Interference. 



b. Performance of the Optimum Receiver 

The above result for the probability of error of the suboptimum receiver 
can be compared to the performance of the optimum receiver using Equation 4.5 and 
Equation 4.6. For the FSK modulation, the optimum receiver performance is obtained 
by first evaluating 



E'(l • p') 
■) 



_Lr f 

“^0 ^ 1 + q-(t);^ 



T 

dt + J Sj-(t) dt ]. 
cT 
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If we let b = 0 and c = 1 as a special case, then Q^lO) = a" 3. Thus we can simplify 
the above equation to obtain 



E'(l - p ) 
-> 



A-T 



2NV 





SNR 



1 + JSR-SNR 
1 + 2JSR-SNR 



)■ 



(5.56) 



When JSR = 0. q(t) = 0. thus 
E'(l-p') SNR 



so that the actual performance being obtained from 



erfc ( y SNR 4 ) 
erfc ( VsNR p 4 ) 

where 

I + JSR-SNR 

P = . 

1 + 2JSR-SNR 

A plot of Pg o ^ function of SNR and JSR is shown in Figure 5.16 with parameters 
b set to zero, and c set to one. 

3. Comparison 

We observe here identical results on suboptimality for FSK as for PRK. That 
is, since a and P are similar in form for both modulation schemes, we find here that 
FSK modulation is no less susceptible to jamming than PSK modulation when using a 
suboptimum receiver to process the signal, the noise, and the jantming. 




JSR = 0 



JSR ^ 0 



(5.58) 
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PROBABILITY OF ERROR 
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Fisure 5.16 Performance of the Optimum Receiver for FSK Modulation 
with Pulsed Wave Interference. 
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VI. CONCLUSIONS 



In this thesis general results have been obtained for the performance of optimum 
and suboptimum (conditional) digital communication receivers operating in the 
presence of (stationary) AWGN and nonstationary AWGN generated via a specific 
model. The general results show that the suboptimum receivers have performance (i.e., 
probability of error, ^ ) that is always inferior to that of their optimum counterparts 
(expressed also as probability of error. ^ ). 

I he many e.xamples worked out in Chapter V for PRK and FSK modulation for 
different nonstationary AWGN interference clearly demonstrate the level of the 
suboptimality of conditional receivers. Perhaps more importantly, the examples 
demonstrate the high degree of vulnerability of conventional receivers to jamming, 
having characteristics similar to the nonstationary AWGN interference a model used. 
The Pg j plots demonstrate that without powerful jamming, that is for relatively low 
JSR values, the receivers can be rendered almost completely ineffective as their error 
probabilities are significantly higher than the lO'"^ BEP military standard. 

The optimum receivers proposed however perform significantly better when 
compared to the suboptimum receivers and in fact in all cases, with sufTicient SNR. are 
able to overcome the effect of jamming. It is important to note however that these 
receivers, in order to be optimum most have knowledge about the nonstationary 
(jamming) interference that would normally not be available. Consequently the actual 
performance of these receivers may in practice be work worse than predicted and 
perhaps worse than the suboptimum receivers. 
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APPENDIX A 

DERIV ATION OF THE APPROXIMATION TO EQ. 5.10 

From Equation 5.10 we must evaluate 



E'(l - p') _ SNR(1 - nr) 
- V3JSR-SXR 



[ tan* V3JSR-SNR + | 



y3JSR-SNR 

cos(47ifTxV3JSR-SNR) 



dx 



If we assume that y3JSR*SNR > > 1, then 



E'(l - p') 
■) 



SNR(l-m‘) , , 

: ■ ■■- - [ tan ** V3JSR-SNR 

v^3JSR-SNR 




g v^3JSR'SNR 



) 



as the above integral can be evaluated in closed form when VJSR’SNR > > 1. 

If we further assume that JSR'SNR is so large that the exponential is nearly 

unity and tan'* v^3JSR'SNR * 7t 2, then 

E'(l-p') _SNR(l-nr);t 
^ 73JSR-SNR 
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APPENDIX B 

DERIV ATION OF THE APPROXIMATION TO EQ. 5.40 



Recall from the discussion that for the suboptimum receiver for FSK modulated 
signal, placing f^ at f^ 2 in the sinusoidal wave interference frequency is most damaging 
to the performance of the receiver. For this reason, f, is set to f^ 2 and Equation 5.40 is 
approximated under this constraint. VV'e obtain 

E'(l - P') _ A^Tn 1 "^^^(1 +COS X) sin- (f^x 2g ^ 

2 n(a" 3) 2;t 1 + n cos x 



where r| was defined before. 

Observe that typically f^ 2f^ < < 1. (For fj = 3.002 MHz and f^ = 3 MHz. f^ 2f^ 
— 1.67 10“^ ). Since (1 + cos x);(l + rj cos x ) is periodic in 2;t, and sin" (fj 2f^) is 

essentially constant for 0 ^ x ^ 27t, we obtain 



Ed - p ) _ A-Tn V 1 f 



1 + COS X ■> 

dx sin" ( — ^ 2;tk ) 

n(a- 3) ' 27t •' Q 1 + q cos x 2f, 



A-Tq 
n{a“, 3) 



I + P ^ ^ 

(1 - p) 5— I X s>n“ ( 



1 - p 



p-l 



2f 



and p is given by Equation 5.30. The finite sum can be evaluated in closed form and we 
finally obtain 



E'(l - p') 
1 



(A-T 2)q 1 + p- [cos (n+ l)(7tfj p ] sin (nTtf^f^) 

V, ->M*‘ . I 

(a- a) • - P" n sin (ttlj ip 



1 

JSR 



2q + 2q7l - q- 



2 + q - q: + (2 + q)yi - q2 



sin (nTtf. f ) 

]{1 - [ cos (n+ l)(nf.,f ) 1 , '■}• 



O 
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For n large, the ternt 



sin (nTtf .r ) 

[cos(n+l)(7Tfj gi ^ 0. 

nTTlc 1^ 

thus we obtain the approximation 

E'(l-p') SNR 

== ( )' 

2 2 

^ 2(1 + JSR-SNR) + 2 /l + 2JSR-SNR 

2 + 5JSR-SNR + 2(JSR-SNRr + (2 + 3JSR-SNR)yi + 2JSR-SNR 
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APPENDIX C 

DERIVATION OF EQUATION 5.44 



For PRK modulation with pulsed wave interference, from Equation 5.44 we 
evaluate 



cT 

J Q(t) Sj^(t) dt = J q^(t) [ 4A^(1 - m') cos^ 2rtf^t ] dt 
T 0 



2a V P-l {E+‘''OT 

= 2A'(l-nr)( +2abc + 2b-) V J' 

^ E^O CT 

P 



'' ■> 
a"c“ 



P -1 



(1 + COS 47tf^t) dt 
sin 47tf^t (C+' 2 )T 



= 4A*(1 - m“)( f-abc + b") T [t h )' ' P 

3 t=o 'Tp 



2 



P-i 



= 2A“(1 - nr)(-^ — +abc + b“) Y (T + [sin 47tf^(C+ '/ 2 )T - sin 4rtf^CT J 2rtg. 



C = 0 



Thus, for the sum term we have 



P -1 

V [ sin 4Ttf(C+ '' 2 )T -sin4rtfCT ] 

C = 0 ' 

p-1 p-1 

= sin2rtfT V cosdttfCT - (1 - cos 2Jtf T ) 'V sin47tf£T 

r p — r p ' r p'^ — r p 

C =0 £=0 



= sin 27tfT cos 2n:f(p-l)T sin 27rf pT , sin 27tf T 

r p ^ p p r p 

- 2sin“ Ttf T sin 2nf (p-l)T sin 27tf pT sin 27tf T 

r p r'^ ' p p r p 



since 



k 

y cos 27tfCx = cos k27tfx sin 27tffk-l)x sin 27tfx 
C = 0 
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and 



k 

^ sin 27rfXx = sin k27rfx sin 2;rf(k-r)x sin 27tfx. 

E = 0 

Thus the sum term 
P 

y [ sin4;tf(£+ iT)T -sin47tf£T ] 

E = 0 

= cos (p-l)2Ttf Tp sin 27tf.Tp - tan ttf^Tp sin (p-l)27tf^Tp sin 2trf.pTp 
= { cos 27tf'pT cos 27tf T + sin 27tf pT sin 27tf T ) sin 27tFpT 

' p r p r P r P P 

- tan ttf T sin 27rf pT (sin 2jtf pT cos 2;rf T - cos 2Ttf pT sin 2;tf T ) 

r p p' P P P r p- 

= (sin 47tf^pTp cos 27tf.Tp) 2 + sin" 2trf^pTp sin 2ttf.Tp 

- sin" 2ttf pT cos 27rfT tan trf'T + (sin 47rf pT sin 2Jtf T tan itf T ) 2 

r^p rp rp' r’^P rp rp' 

= (sin 47if cT cos 27Tf cT p) 2 + sin" 27tf cT sin 2trf cT p 

- sin" 2trf'j.cT cos 27tf^cT p tan Jtfj.cT p + (sin 4trf^cT sin 27tf cT p tan ;rrj.cT p) 2 
where pTp = cT. Since sin 2ttf cT/p tan ttf cT'p = 2sin" ttf cT p, 

P '1 

y [ sin 47tf(e+ !'i)T - sin 47tf £T ] 

r' ^p rp-* 

£=0 

= [ sin dtif^cT (cos 2ttf cT p + 2sin" 27tf cT p) ]'2 

+ sin" 27rf cT (sin 27tf^cT;p - cos 2Ttf^cT p tan rtf^cT p). 
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Observe that the terms in tlie lirst parenthesis is unity and in the second parenthesis is 
tan Ttf^cT p. thus 



P -1 

^ [ sin 4;tr(£+ ‘^ 2 )! - sin 47tf;£T ] 
£ = 0 



= sin" 27tf^T tan ttfcT p + (sin 4rtfcT) 2. 



The remaining sum term is 



P -I 

y (T +1 27tf) = cT(l + l'2rtf;cT). 
£ = 0 



Thus from above results we finally obtain 



j q"{t)s^"(t)dt = 2A"T( 1-nr) — (a"c"+ 3abc+ 3b")* 

T ^ 



[ 1 H (sin 4nf cT + 2sin" 27tf cT tan ttf cT p ) ]. 

^ 4rtfrcT ' ^ r H / J 
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APPENDIX D 

DERIVATION OF EQUATION 5.50 

The sum term in Equation 5.50, namely 
P -1 

y [ sin Jtkc(C+ !' 2 )-P - sin TtkcC p ] 

£ = 0 

yields 

P -* 

= [ sin TtkcC p cos ttkc 2p + cos Ttkcf p sin ttkc 2p - sin TtkcC p j 

£ = 0 

p-1 p-1 

= sin 7tkc;2p T cos ttkcC p - (1 - cos Ttkc 2p) T sin ttkcE p 

£=o 

= [ sin JTkc, 2p cos (p-l )7tkc,'2p sin rtkcp 2p j sin ttkc, 2p 
- [ 2sin‘ Ttkc 4p sin (p-l)7tkc/2p sin Ttkcp 2p j. sin Ttkc'2p 

= cos (p-l)7tkc'2p sin Ttkc'2 - tan Ttkc,4p sin (p-1) Ttkc 2p sin Ttkc 2 

= sin Ttkc, 2 [cos (p-l)Ttkc'2p - tan Ttkc,4p sin (p-l)Ttkc,2p ] 

= sin Ttkc 2 [ sin Ttkc 2 (sin Ttkc, 2p - tan Ttkc,'4p cos Ttkc'2p) 

-I- cos Ttkc/2 (cos Ttkc,'2p -4 tan Ttkc;4p sin Ttkc,'2p) ] 

= sin Ttkc, '2 (sin Ttkc. 2 tan Ttkc'4p -i- cos Ttkc 2) 

= sin* Ttkc 2 tan Ttkc, 4p -I- 2sin Ttkc. 



73 






LIST OF REFERENCES 



1. Van Trees. H. L. and Wilev, J., Deieciion, Estimaiion, and Modulation Theorv. 
Chapter 4, 196S. 



Whalen, A. D., Detection of Signals in .\'oise. Chapter 6, Academic Press, 1971. 



Penfield, P.. Jr., "Fourier Coefilcients for Power-Law Devices," Journal of the 
Franklin Institute, pp. 107-122, February’ 1962. 



74 



INITIAL DISTRIBUTION LIST 



No. Copies 



I. Defense Technical Information Center 2 

Cameron Station 
Alexandria. Vireinia 22304-6145 



2. Librarv. Code 0142 2 

Naval Postgraduate School 
Monterey, California 93943-5002 



3. Department Chairman. Code 62 

Department of Electrical and Computer Engineering 
Naval Postgraduate School 
Monterey, California 93943 



4. Professor D. C. Bukofzer. Code 62Bh 

Department of Electrical and Computer Engineering 
Naval Postgraduate School 
-Monterey, California 93943 



5. Professor G. A. Mvers. Code 62Mv 2 

Department of Electrical and Computer Engineering 
Naval Postgraduate School 
Monterey, California 93943 



6. Kim. Young Joo 5 

213- 901. Han- sin gong yung 6th Apartment 
Ban po-Dong, Gang nam-Gu 
Seoul. 130-09, Korea 



7. Hwang. Chun Sik 1 

Departmemt of Electrical Engineering, 

Militarv Academv, P.O.Box 77, Gong neung-Dong 
Do bohg-Gu. Seoul, 130-09, Korea 



75 






-s 7 ^ 

r 



i'AiJX'rr.tDW.jri; Y , UALlFORi.'] 



iThesis 

>K42125 




22105: 



Kim 

Performance of digital 
communication receivers 
operating in the presence 
of white Gaussian nonsta- 
tionary noise. 



2?i051 



Thesis 
K42125 
c. 1 



Performance of digital 
communication receivers 
operating in the presence 
of white Gaussian nonsta- 
tionary noise. 



